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Abstract 

A nunnber of recent studies have shown that intraspecific genetic variation of 
plants nnay have a profound effect on the herbivorous connnnunities which 
depend on thenn. However less is known about the relative innportance of 
intraspecific variation connpared to other ecological factors, for exannple 
environnnental variation or the effects of herbivore dannage. We randonnly 
selected 22 Betula pendula genotypes fronn a local population (< 0.9 ha), 
cloned thenn and planted cloned seedlings on two study sites separated at a 
regional scale (distance between sites about 30 knn) to exannine an insect 
connnnunity of 23-27 species on these genotypes. B. pendula genotypes did not 
differ in their species richness, but the total nnean abundance and the structure 
of the insect herbivore connnnunity was significantly affected by the genotype, 
which could account for up to 27% of the total variation in connnnunity structure. 
B. pendula genotype accounted for two to four tinnes nnore variation in the 
arthropod connnnunity structure than did environnnental (block) variation on a 
local scale, while on a regional scale, genotypic and environnnental (site) 
variation accounted for 4-14% of the arthropod connnnunity structure. The 
genetic effects were nnodified by environnnental variation on both a local and 
regional scale over one study year, and locally, the largest part of the variation 
(38%) could be explained by the genotype x environnnent (block) interactions. 
Suppression of insect herbivores during one growing season led to changed 
arthropod connnnunity structure in the following growing season, but this effect 
was nnininnal and could explain only 4% of the total variation in insect 
connnnunity structure. Our results suggest that both genetic and environnnental 
factors are innportant deternninants of the connnnunity structure of herbivorous 
insects. Together these nnechanisnns appear to nnaintain the high diversity of 
insects in B. pendula forest ecosystenns. 



Article Status Summary 
Referee Responses 



Referees i 


1 


2 


v1 


□ 


u 


published ! 


report 


report 


31 Jan 2014 1 






1 Patrick Tobin, United States Department 
of Agriculture USA 



2 Judith l\/lyers. University of of British 
Columbia Canada 

Latest Comments 

No Comments Yet 



Page 1 of 13 



FlOOOResearch 



FIOOOResearch 2014, 3:34 Last updated: 24 FEB 2014 



Corresponding author: Tarja Silfver (Tarja.h.silfver@gnnail.conn) 

How to cite tiiis article: Silfver T, Rousi M, Oksanen E etal. (2014) Genetic and environmental determinants of insect herbivore community 
structure in a Betula pendula population [v1 ; ref status: indexed, http://f1 000r.es/2pd] FIOOOResearch 2014, 3:34 (doi: 
1 0.1 2688/f 1 000research.3-34.v1 ) 

Copyright: © 201 4 Silfver T et al. This is an open access article distributed under the terms of the Creative Commons Attribution Licence, which 
permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The author(s) is/are employees 
of the US NIH and therefore any publishing licenses are also subject to the terms of the NIH Publishing Agreement and Manuscript Cover Sheet. 
Data associated with the article are available under the terms of the Creative Commons Zero "No rights reserved" data waiver (CCO 1 .0 Public 
domain dedication). 

Grant information: This study was funded by the Academy of Finland (project 78743, granted to EO), and funding granted to TS by the Graduate 

School in Forest Sciences, Finnish Concordia Fund and Kone Foundation (grant 2-712) is also acknowledged. 

The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. 

Competing Interests: No competing interests were disclosed. 

First Published: 31 Jan 2014, 3:34 (doi: 10.1 2688/f 1000research.3-34.v1) 
First Indexed: 24 Feb 2014, 3:34 (doi: 10.1 2688/f 1000research.3-34.v1) 



Page 2 of 13 



FIOOOResearch 2014, 3:34 Last updated: 24 FEB 2014 



Introduction 

Genetic variation within one species can affect the structure and 
dynamics of associated communities and entire ecosystems^ -. This 
may be considerable, especially for keystone species, such as forest 
trees, which serve as food and habitat for numerous primary con- 
sumers. A vast number of studies have already shown that arthro- 
pod communities respond to genetic differences among individual 
plants within interspecific hybridizing complexes (e.g. Eucalyptus^\ 
Salix^, Populus^, Quercus^') or specific genotypes within species 
(e.g. Oenothera biennis\ Eucalyptus globulus^, Solidago altissima^, 
Populus angustifolia^^^). However, it has recently been argued that 
the role of plant genetic variation in structuring arthropod com- 
munities has been considerably inflated due to the common meth- 
odological flaw that genotypes are collected from diverse and often 
distant environments, which maximizes genetic variation, whilst 
experiments are performed in a single common garden where envi- 
ronmental variation is minimized^ Indeed, when this mismatch 
in scale was avoided in the experimental design, spatial processes 
relegated host plant genotype to a secondary role in structuring 
insect communities of Quercus robur L.^\ Whether this applies to 
all systems is, however, not yet known. 

Genes encounter a range of environments in nature and it has long 
been recognized that genetic determination of plant susceptibility to 
a herbivorous insect depends on environmental context^^. However, 
most studies that have examined the role of genotype x environ- 
ment interactions in the abundance and distribution of herbivorous 
species, have used only one or a few closely related herbivore spe- 
cies (e.g.^''"^^), and much fewer studies have examined genotype x 
environment interactions in a community context^' ^^'^"^ -^^ It is well 
recognized that we know too little of the relative importance of 
intraspecific genetic variation compared to other ecological factors 
that also influence multi-trophic communities and ecosystem 
processes^ ^ Thus, the examination of genotype x environment inter- 
actions in a community context may be essential for improving our 
knowledge in the developing field of community genetics. 

Silver birch (Betula pendula Roth) is an ideal tree species in which 
to examine the mechanisms of plant-herbivore interactions and the 
community-level consequences of trait variation, because the species 
shows remarkable genetic variation in its resistance to herbivores^ 
In addition, the genetic variation of secondary metabolites^^', nutri- 
ent concentrations^^, and phenological traits-^ -'^ of B. pendula are 
known to be substantial, and all these traits are known to affect 
herbivores and higher trophic level interactions^'^^^'^^ Most of the 
studies that have been conducted using B. pendula have used 
genotypes that were originally randomly selected from a local 
B. pendula population, i.e. from a naturally regenerated forest stand 
< 0.9 ha. None of these earlier studies have, however, investigated 
the within-population genotypic variation in B. pendula insect her- 
bivore species richness and community composition. We cloned 
22 B. pendula genotypes, planted them in two common gardens sep- 
arated at a regional scale (distance between sites about 30 km), and 
studied the relative importance of genetic variation in community 
patterns, comparing both local and regional environmental varia- 
tion. In addition, we examined how strongly herbivores themselves 
can modify arthropod communities associated with B. pendula by 
suppressing herbivores from half of the saplings over one growing 



season in one common garden and surveying their arthropod com- 
munities the following season. 

Materials and methods 

Plant material and study sites 

The 22 different genotypes of B. pendula were cloned during spring 
1998 from randomly selected B. pendula trees taken from a natu- 
rally regenerated B. pendula - B. pubescens Ehr forest in Punkaharju, 
southeastern Finland (61°48' N, 29° 18' E), to study genetic varia- 
tion in phenology, growth, reproduction and resistance-related traits 
among individual birch trees-^ Sampling was stratified random 
sampling: six spots where forest lift could be transferred were first 
selected around the forest, and 2-5 trees within the reach of forest 
lift in each spot were then randomly (by throwing a coin) selected 
for our study purposes. B. pendula is predominantly a sexual spe- 
cies, but genotypes can be cloned for study purposes or for planta- 
tions using standard tissue-culture methods^^. Cloned B. pendula 
saplings were planted at the growing sites (i.e. common gardens, 
each approximately 0.25 ha) in June 1999 to find out the degree 
to which the genotype and environment affect birch traits and to 
test how genotypes differ in their response to the environment-^'. 
The Kuikanniitty study site (61°47' N, 29°2V E) is an abandoned 
cultivated field and the Parikkala study site (61°36' N, 29°36' E) is 
Myrtillus type forest^l Soil type was defined as fine sandy till for 
both sites-^\ The distance between these sites was around 30 km and 
they were situated at approximately the same altitude (Kuikanniitty 
79 m and Parikkala 93 m above sea level). Thus, the mean sum- 
mer (June- August) temperatures were very similar at these sites: 
in 2002 mean temperatures were 17.6°C and 17.9°C and in 2003 
they were 15.9°C and 15.6°C in Kuikanniitty and Parikkala, respec- 
tively. Both study sites were divided into six blocks, each of which 
included four saplings from each genotype. To prevent edge effects, 
the experimental saplings were surrounded by one row of extra sap- 
lings. From each block, one of the four saplings of a total of 22 
genotypes was randomly selected for the present study in order to 
have six replicates per genotype. 

In addition, we collected additional data from Kuikanniitty in 2003 
to investigate the effect of previous insect herbivory on insect com- 
munity structure and abundance, and surveyed one extra sapling 
from each block and genotype. These extra saplings were protected 
from insect herbivory in the previous growing season by regular 
sprayings with synthetic pyrethrin-\ which has no direct or side 
effects on the growth or chemistry of birch seedlings^^. 

Measuring insect abundance and species richness 

The insect herbivore community of each sapling was assessed 
by surveying the abundance of 23 (in Parikkala 2002) or 27 (in 
Kuikanniitty 2002-2003, and Parikkala 2003) insect taxa from diverse 
orders (Lepidoptera, Hymenoptera, Coleoptera, Diptera, Hemiptera; 
Table 1). These taxa were generally the most abundant taxa in both 
sites. However, species that were rare in both sites were included 
in the surveys as well. Species identifications were undertaken fol- 
lowing Saalas^^ species identification guide, using several web pages 
(http://www.funet.fi/pub/sci/bio/life/insecta/index.html; http://www. 
leafmines.co.uk/index.htm; http://www.bladmineerders.nl/; http:// 
www.nrm.se/) with the assistance of specialists. Euceraphis betulae 
eggs were counted from the side of twelve (2002) or eight (2003) 
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Table 1. Description of the 27 taxa surveyed for their abundance among 22 genotypes in Kuikanniitty and Parikkala field 
experiments 2002 and 2003. 







2002 




2003 








Kuikanniitty 


Parikkala 


Kuikanniitty Parikkala 


Taxa ^unnmm^^^^^^ 




1 Total number of insects/damage counted 




Lepidopteran miners/rollers 


Gracillaridae 1 (miner) 


Phyllonorycter cave 1 la 


282 


123 


53 


34 


Gracillaridae 2 (miner) 


Phyllonoryctersp. 1 


61 


26 


12 


4 


Gracillaridae 3 (miner) 


Phyllonorycter sp. 2 


19 


7 


3 


3 


Gracillaridae 4 (miner) 


Parornix betulae 


114 


42 


40 


10 


Gracillaridae 5 (miner) 


Parornixsp. 


30 


11 


20 


0 


Eriocranidae (miner) 


Eriocrania sp. 


536 


2007 


746 


2374 


Pyralidae (roller or tier) 


tentatively Euzophora fuliginosella 


67 


77 


135 


142 


Tortricidae (galler) 


Epinotia tetraquetrana'' 


159 


136 


159 


136 


Nepticulidae (miner) 


Stigmella sp. 1 


40 


53 


7 


1 


Incurvanidae (miner) 


Phylloporia bistrigella 


125 


6 


30 


8 


Geometridae (roller or tier) 


Rheumaptera hastata 


11 


6 


4 


0 


Gelechiidae (roller or tier) 


tentatively Teleiodes sp. 


87 




211 


37 


Mircolepidoptera 1 (roller or tier) 




64 


65 


188 


60 


Lepidoptera 1 (roller or tier) 




8 


2 


3 


0 


Lepidoptera 2 (miner) 




12 




7 


3 


Lepidoptera 3 (roller or tier) 




0 


1 


13 


1 


Lepidoptera 4 (miner) 




142 


7 


152 


82 


Coleopterans 


Attelabidae (roller) 


Deporaus betulae 


62 


14 


157 


133 


Curculionidae (miner) 


Orchestes rusci 


54 


127 


23 


12 


Hymenopterans 


Tenthredinidae 1 (miner) 


tentatively Fenusa pumila 


149 


109 


66 


59 


Tenthredinidae 2 (leaf feeder) 


Hemichroa austral is 


167 




52 


18 


Tenthredinidae 3 (leaf feeder) 


Croesus septentrionalis 


108 


7 


34 


0 


Cimbicidae (leaf feeder) 


Trichiosoma sp. 


6 


2 


2 


0 


Dipterans 


Agromyzidae (miner) 1 


Agromyza ainibetulae 


24 


11 


33 


6 


Cecidomyiidae (miner) 1 




0 


0 


1 


19 


Hemipteran 


Aphidoidea (sap sucker) 


Euceraphis betulae 


996 


2640 


40 


114 


Heteropteran 


Heteropteran 1 (sap sucker) 




92 




284 


466 



E. tetraquetrana counts represent the damage during the whole lifetime of the saplings (see Materials and methods). Note also that years are not directly 
comparable because of the changed sampling protocol between years. 
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topmost buds in April before budburst. However, in 2002 E. betulae 
eggs were counted from two saplings per genotype per block (sum 
of the eggs on the sides of 24 buds was used in the analysis), 
because regular sprayings with synthetic pyrethrin on the other sap- 
ling was started only after egg counts in both sites. Trichiosoma sp. 
pupae were counted in April/May when the timing of budburst of 
the same saplings was observed (Possen, submitted manuscript). 
The abundance of Eriocrania sp. was determined at the end of June, 
Deporaus betulae at the beginning of July and Heteropteran 1 (sap 
sucker) in August. Croesus septentrionalis larval colonies and the 
number of larvae in each colony were recorded along with Eriocrania 
and D. betulae measurements in both years. The abundance of all 
other insects were determined indirectly by counting damaged 
leaves at the beginning of September in both years, since the dam- 
age caused by most of the surveyed taxa remained identifiable for a 
long time after the initial damage. 

In general, the insect abundance in 2002 was determined by sur- 
veying the whole sapling. The mean height of these saplings at the 
end of 2002 was 253 ± 4.3 cm (mean ± SE) in Kuikanniitty and 
227 ± 3.8 cm in Parikkala. Because B. pendula genotypes differ in 
their height and diameter growth-^ and large saplings may harbor 
more insects than smaller saplings, we determined the whole sap- 
ling "surface area" and used it as a covariate hereafter called "size 
index" in statistical analysis. Surface area was determined by pho- 
tographing each sapling sideways from their southern side against 
a white background, converting the picture to a black and white sil- 
houette picture in Adobe Photoshop 7.0 and determining the num- 
ber of black pixels (i.e. leaf and branch area) within the picture. The 
number of pixels was converted to m^ using the number of pixels 
of a known area as a reference. The amount of pixels significantly 
(p < 0.001) explained over 73% of the sapling volume [Y = (3.14 
* {base diameter/2}^ * height)/3] in both sites. The abundance of 
Phyllonorycter cavella, Phyllonorycter sp. 1, Parornix betulae and 
Parornix sp. was not examined on the whole sapling, but was 
determined as the damage (i.e. number of mines per each species) 
found within a period of 30 seconds. The period of time (30 sec) 
was chosen so that even the smallest saplings had leaves uncounted 
when the time was up. 

Since the method of assessing herbivore abundance/resistance by 
time counts has been successfully used in the past^^'^^' we decided 
to use time counts to determine the abundance of almost all taxa 
(except E. betulae, Trichiosoma sp. and C. septentrionalis) in 2003. 
The same person undertook all surveys. The abundance of easily 
visible damage (large mines and rolls) of Eriocrania sp., Z). betulae 
and Heteropteran 1 were determined as the number of damaged 
areas found within a period of 30 seconds. Epinotia tetraquetrana 
"knobs" in the branches of the saplings were counted within a 
period of 20 seconds in 2003 starting at the top of the tree. Since 
the "knobs" in the branches remain visible for years and we did 
not separate different year's growth while surveying, the values 
represent the accumulation of E. tetraquetrana damage during the 
last few years. Therefore the same values were used in both years' 
insect community analyses. The abundance of all other 20 taxa in 
2003 was determined within a single time count of each sapling at 
the beginning of September. To examine a similar proportion of each 
sapling, they were divided into three size categories according to their 



height and number of leaves. Small saplings (average height 2.8 and 
3.2 m in Parikkala and Kuikanniitty, respectively) were surveyed 
for 30 seconds, average sized saplings (3.5 and 3.9 m in Parikkala 
and Kuikanniitty, respectively) for 60 seconds and large saplings 
(4.5 m in Kuikanniitty, large saplings were not found in Parikkala) 
for 120 seconds. Surveying time was used as a covariate called size 
index in statistical analysis. 

Data analyses 

All multivariate analyses were performed with Primer 6 (Primer-E 
Ltd, United Kingdom). The full data matrix consists of the abun- 
dance of 23-27 (23 in Parikkala 2002) insect species in 264 saplings 
(22 genotypes, 6 blocks, 2 sites) that were surveyed in two consecu- 
tive years. All surveyed insect species were included in the statisti- 
cal analysis when sites were tested separately, but those four species 
that were not surveyed in Parikkala 2002 were excluded also from 
Kuikanniitty 2002 data when sites were compared. Arthropod com- 
munity composition data was analyzed using non-parametric mul- 
tivariate analysis of variance (PERMANOVA), which is well suited 
to non-normal ecological data such as ours^^'^^. Years were analyzed 
separately in all statistical tests, because of the changed sampling 
protocol between years (surveying the whole tree in 2002, using 
time counts in 2003). All data was fourth root transformed prior 
to analysis to reduce differences between common and rare spe- 
cies. The semimetric Bay-Curtis distance, which generally seems 
to provide the most meaningful measure of dissimilarity in ecologi- 
cal community structure ''^, was used to calculate distances between 
each pair of observations. The resulting distance matrix was used 
to obtain p-values using a random subset of 4999 permutations in 
PERMANOVA. The permutation method was permutation of residu- 
als under a reduced model. The statistical model was designed to 
test the effect of genotype, site, block (nested within site) and the 
interaction of genotype x site using sapling size index (sapling 
surface area in 2002 and surveying time in 2003, see above) as a 
covariate. Site was treated as a fixed factor and block and geno- 
type as random factors in the model. In addition to these analy- 
ses, we separately tested the effect of genotype and block on insect 
assemblages in each site and year to calculate the proportion of 
variance explained by B. pendula genotype and local environment 
(i.e. replicated block). Additional data collected from those saplings 
that were protected from insect herbivory in the previous growing 
season in Kuikanniitty 2003, were combined with the Kuikanniitty 
2003 non-treated sapling data prior to analyzing the effects of insect 
removal, block and genotype, and their two-way interactions with 
the insect assemblages with PERMANOVA. Sapling size index was 
used as a covariate. 

To visualize the multivariate patterns among observations, non- 
metric multidimensional scaling (nMDS) was performed on the 
Bay-Curtis distances. The distance among centroids for groups of 
samples was determined prior to nMDS to increase clarity, e.g. when 
the whole data was visualized we had 88 genotype- site-year points 
(22 genotypes in 2 sites over 2 years) instead of 528 genotype- 
block- site-year points. To visualize the effect of genotype in indi- 
vidual site and year, we separately determined the distance among 
genotype centroids in each site and year and produced one nMDS 
plot from each of these "environments". Additional Kuikanniitty 
2003 data combined with Kuikanniitty 2003 raw data was used to 
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visualize the effect of insect removal on insect assemblages using 
nMDS on the genotype centroids of those saplings that were either 
protected from herbivory or grown under natural herbivory. 

Species richness (number of species/sapling) and total mean abun- 
dance (number of herbivores/sapling) was statistically tested by anal- 
ysis of covariance using SPSS 20.0.0.1 (IBM SPSS Statistics) General 
Linear Models (GLM) procedure. Those four species that were not 
surveyed in Parikkala 2002 were excluded also from Kuikanniitty 
species richness and total mean abundance calculations to better 
enable site comparisons. Genotype and block (nested within site) 
were treated as random factors and site as a fixed factor in the sta- 
tistical model while sapling size index was used as a covariate. 
Additional Kuikanniitty 2003 data combined with Kuikanniitty 2003 
basic data was used to analyze the effects of insect removal, block 
and genotype, and their interactions with the species richness and 
total mean abundance. Genotype and block were treated as ran- 
dom factors and insect removal as a fixed factor while sapling size 
index was used as a covariate. Total mean abundance was log(x+l)- 
transformed to equalize the error variances across groups in both 
analyses. 

Results 

Study years and sites were distinctly grouped apart into two- 
dimensional ordination space, when the genotype centroids of different 
years and sites were analyzed using nMDS (Figure 1). The MANOVA 
Table in tum, shows that sites had statistically significantly different 
insect species community composition in both years (Table 2). Sites 
were also clearly different in their total mean abundance (Table 3) 
and species richness (p < 0.008 for the site effect in species rich- 
ness): the forest site of Parikkala had a 49-78% higher total mean 
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Figure 1. Non-metric UDS plot of insect assemblages of 23 
species colonizing 22 B. pendula genotypes in Kuikanniitty and 
Parikkala 2002 and 2003. Stress is 0.14, which indicates a good 
representation of the data in two-dimensional ordination plot. Each 
point is a centroid of six replicates. Numbers in the centre of the 
markers are genotype identification numbers. 



abundance, but 18-25% lower species richness than the abandoned 
field site of Kuikanniitty in 2002-2003, respectively. These findings 
indicate that each year and site had significantly different herbivo- 
rous insect assemblages, thus creating different biotic environments. 

Genotypic variation and genotype x environment interactions 

B. pendula genotypes were significantly different in their insect 
species community composition in both study years (Table 2). In 
2002, regional scale environmental (site) variation explained more 
of the total variation in species composition than the genotype (13.9 
and 8.0%, respectively), while in 2003 the genotype explained more 
of the total variation than the site (12.1 and 3.8%, respectively). 
Significant genotype x site interaction, which explained 8.2% of 
the total variation, was found only in 2003. When the sites were 
tested separately in both years we found that the effect of geno- 
type was significant in Kuikanniitty 2002 and both study sites in 
2003 (Table 4, Figure 2). B. pendula genotype could account for 
15.8-27.0% of the total variation in community structure, while 
local scale environmental (block) variation explained 5.9-7.6% of 
the total variation in community structure (Table 4, Figure 2). 

B. pendula genotypes also significantly differed in their total mean 
abundance of herbivores (mean number of herbivores/sapling): the 
total mean abundance of the most susceptible genotype was 5.4- 
and 3.2-fold compared to the total mean abundance of the most 
resistant genotype in Kuikanniitty and Parikkala 2002, respectively 
(Table 3, Figure 3). In 2003, only the genotype x site interaction 
was statistically significant, which indicates that the genotype 
effect strongly depended on the study site. Indeed, when we tested 
the study sites separately, genotype effect was significant only in 
Parikkala (ANCOVA: Parikkala F^^ ,^^='^•'2.9, p=0.003; Kuikanniitty 
^21 io3~^'^^' P=0-103). The species richness (number of insect spe- 
cies/sapling) was not significantly affected by the B. pendula geno- 
type or genotype x site interactions in either year (p>0.134). 

Local scale genotype x environment interaction (i.e. the interaction 
of genotype x replicated block) was studied in Kuikanniitty 2003. 
Insect species community composition was significantly affected by 
both genotype and genotype x block interaction (Table 5). Genotype 
variation explained 10.6% and genotype x block variation 38.0% 
of the total variation in insect community composition, indicating 
that genotype effect is also strongly affected by local scale environ- 
mental variation. Total mean abundance or species richness was not 
affected by genotype or genotype x block interaction (p>0.097). 

Effects of the previous year's herbivory on insect communities 

Previous year herbivory changed the insect community composi- 
tion of B. pendula saplings (Table 5). The genotype centroids of 
those saplings that were either subjected to natural herbivory or 
protected from it were located on the opposite sides of the two- 
dimensional nMDS ordination plot, although overlapping is evi- 
dent (Figure 4). Previous year herbivory did, however, explain 
only 4.4% of the total variation in insect community composition. 
Total mean abundance was affected by the previous year's herbivo- 
ry as well, but species richness was not (ANCOVA: effects of insect 
removal on total mean abundance F^ ^ ^8-^4.6, p=0.002 and species 
richness p>0.829). 
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Table 2. Non-parametric MANOVA table of the effects of genotype and site on insect 
herbivore community structure on B. pendula saplings in 2002-2003. Sapling size 
index, wliicli is a measure of lieiglit and number of leaves (see material and methods), was 
used as a covariate. 



Insect herbivore community 2002 Insect herbivore community 2003 





df 


SS 


F 


P 


df 


SS 


F 


P 


Genotype 


21 


27700 


1.26 


0.04 


21 


33294 


2.04 


< 0.001 


Site 


1 


47798 


16.99 


< 0.001 


1 


10319 


7.62 


< 0.001 


Block (Site) 


10 


17826 


1.70 


0.002 


10 


18013 


2.31 


< 0.001 


G X S 


21 


22242 


1.01 


0.466 


21 


22472 


1.37 


0.006 


Size index 


1 


4890 


4.65 


< 0.001 


1 


8084 


10.38 


< 0.001 


Residual 


208 


218510 






208 


162000 






Total 


262 


344260 






262 


274880 







Table 3. The ANCOVA table of the effects of genotype, block and site on 
total mean abundance of herbivores (log[x+1] transformed) of B. pendula 
saplings in 2002-2003. Sapling size index, which is a measure of height and 
number of leaves (see material and methods), was used as a covariate. 

Mean abundance 2002 Mean abundance 2003 | 





df 


SS 


F 


P 


df 


SS 


F 


P 


Genotype 


21 


1.10 


2.82 


0.011 


21 


0.26 


1.08 


0.435 


Error 


21 


0.39 






21 


0.24 






Site 


1 


1.13 


14.6 


0.004 


1 


0.59 


46.7 


< 0.001 


Error 


9.6 


0.75 






26.7 


0.34 






G X S 


21 


0.39 


1.00 


0.470 


21 


0.24 


1.81 


0.019 


Error 


208 


3.88 






208 


1.30 






Block (Site) 


10 


0.77 


4.11 


< 0.001 


10 


0.16 


2.55 


0.006 


Error 


208 


3.88 






208 


1.30 






Size index 


1 


0.30 


16.3 


< 0.001 


1 


0.23 


36.4 


< 0.001 


Error 


208 


3.88 






208 


1.30 







Table 4. Non-parametric MANOVA table of the effects of genotype and block on insect herbivore 
community structure on 8. pendula saplings in Kuikanniitty and Parikkala 2002-2003. Sapling 
size index, which is a measure of height and number of leaves (see material and methods), was 
used as a covariate. 





Insect herbivore community 2002 


Insect herbivore community 2003 




df 


SS 


F 


P 


df 


SS 


F 


P 


Kuikanniitty 


Genotype 


21 


33778 


1.29 


0.018 


21 


25678 


1.33 


0.013 


Block 


5 


10803 


1.73 


0.006 


5 


10305 


2.25 


< 0.001 


Size index 


1 


4468 


3.58 


0.001 


1 


5948 


6.48 


< 0.001 


Residual 


103 


128690 






103 


94495 






Total 


130 


181760 






130 


138930 






Parii<i<ala 


Genotype 


21 


16771 


1.08 


0.297 


21 


29083 


2.16 


< 0.001 


Block 


5 


8067 


2.17 


< 0.001 


5 


7754 


2.41 


< 0.001 


Size index 


1 


1735 


2.34 


0.033 


1 


2811 


4.38 


0.001 


Residual 


104 


77278 






104 


66834 






Total 


131 


105910 






131 


107680 
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Figure 2. Non-metric MDS plot of insect assemblages of 23 
(Parikkala 2002) or 27 species colonizing B. pendula genotypes in 
a) Kuikanniitty 2002, b) Parikkala 2002, c) Kuikanniitty 2003 and d) 
Parikkala 2003. Each point is a centroid of six replicates. Numbers 
in the centre of the markers are genotype identification numbers. 
White circles denote genotypes in Kuikanniitty, black circles denote 
genotypes in Parikkala. Stress values >0.2 indicate that this data 
may be better visualized with more dimensions (stress for three- 
dimensional solutions varied between 0.14 to 0.16). 

Correlations between species 

The associations between insect species across genotypes in dif- 
ferent sites and years seemed to be based on random associations, 
since we found only one correlation that was significant after 
sequential Bonferroni correction"^". An unidentified gallery mine 
(Lepidoptera 4) and E. fuliginosella were correlated across geno- 
types in Kuikanniitty 2002 (Pearson's correlation; in 2002, r = 0.88, 
n = 22, p < 0.0001; in 2003, r = 0.545, n = 22, p = 0.009), but not in 
Parikkala (p>0.199). 



Community structure of insect lierbivores on different genotypes 
of silver bircli {Betula pendula) 

3 Data Files 

http://dx.d0i.0rg/l 0.6084/m9.f igshare.91 5332 

V ^ 

Discussion 

Our results provide evidence that genetic variation within a natural 
B. pendula population can modify the structure of the arthropod com- 
munity even though all genotypes supported similar insect species 
richness. Genetic variation in phenotypic plasticity, however, seemed 
to be the major factor affecting the abundance and structure of the 
insect herbivores associated with this tree species, because geno- 
type effect was often dependent on the environmental variation at 
both regional (Table 2 and Table 3) and local scales (Table 5). Those 
B. pendula genotypes that were used in our study should give unbi- 
ased estimates of the true variance that is present in B. pendula popu- 
lations, since we chose them randomly from one naturally regenerated 



population stand (< 0.9 ha) in eastern Finland, where this Eurasian 
deciduous tree species is particularly abundant^^. By contrast, we 
might have exaggerated the role of regional environmental varia- 
tion and genotype x environment (site) interactions by planting our 
genotypes on two rather different areas (open forest and abandoned 
field, areas that are typically rapidly colonized by B. pendula) at a 
much larger scale (70,000 ha). Therefore, it is not surprising that the 
importance of the genetic variation in structuring insect herbivore 
communities of B. pendula decreased from 15.8-27.0% (of variation 
explained) to 8.0-12.1% with increasing spatial scale in our study. 
Other studies have also found that while the effect of a genotype can 
be clear on local scales (within common gardens), it may be partially 
swamped by environmental variation on larger scales^'^\ 

It has been argued that, because host plant genotypes have often been 
collected from large geographic areas and studied within the con- 
fines of a single common garden, the role of the host plant genotype 
in arthropod community patterns has been largely overestimated^^. 
Indeed, Tack et aV^ showed that spatial processes dominated 
genetic effects when genotypes of Q. robur were collected at the 
same local (500 ha) or regional (1 million ha) scale as that where 
experiments were conducted, and thus, in real landscapes, spatial 
impacts might relegate host plant genotype to a minor role. Our 
results, however, suggest otherwise, because genotype explained 
about three times more of the total variation in insect herbivore 
community structure than local environment (block) in both sites 
(Table 4), and the scale of our common garden(s) was approxi- 
mately the same as the scale of that where genotypes were collected 
(< 0.9 ha). In addition, on a regional scale, genetic and environ- 
mental effects explained similar proportions of the total variation in 
arthropod community structure (Table 2), even though we might 
have inflated the role of the environment in our study. This dis- 
crepancy in our results might perhaps be attributed to the differ- 
ence in the distribution of these wind-pollinated tree species: the 
populations of Q. robur are strongly fragmented and grow at the 
northern margin of the species' European distribution in southern 
Finland (where Tack et al}^ conducted their experiments), while 



□ Kuikanniitty ■Paril<kala 




3 23 5 2 18 14 19 7 20 4 26 30 6 16 25 15 8 22 9 24 12 17 



Genotype 

Figure 3. l\/lean abundance of insect herbivores (±SE) among 
B. pendula genotypes in Kuikanniitty and Parikkala study sites 
in 2002-2003. White bars: Kuikanniitty, black bars: Parikkala. 
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Table 5. Non-parametric MANOVA table of the effects of 
genotype, block and previous year Insect removal on Insect 
herbivore community structure among B. pendula saplings 
In Kulkannlltty 2003. Sapling size index, which is a measure of 
height and number of leaves (see material and methods), was 
used as a covariate. 





df 


SS 


F 


P 


Genotype 


21 


30358 


1.39 


0.004 


Insect removal 


1 


12525 


7.11 


< 0.001 


Block 


5 


10930 


2.13 


< 0.001 


G X IR 


21 


15557 


0.90 


0.743 


G X B 


104 


109030 


1.27 


0.002 


IR X B 


5 


6031 


1.46 


0.069 


Size index 


1 


2078 


2.52 


0.023 


Residual 


103 


84964 






Total 


261 


286640 







B. pendula has a wider and more continuous distribution over the 
whole of Finland, apart from Lapland. Q. robur populations exhibit 
higher geographic differentiation estimates, F ^ 0.032 for B. pendula 
and 0.066 for Q. robur^^ '^^, which means that the gene flow among 
B. pendula populations is two times higher than among Q. robur 
populations, and thus local B. pendula populations might express 
a larger amount of genetic variation than populations of Q. robur. 

We found that insect herbivore communities can be affected by 
both local and regional genotype x environment interactions, at 
least in some years. But why do B. pendula genotypes support dif- 
ferent insect communities in different environments? It is possible 
that resistance traits of the genotypes are changed due to differ- 
ences in abiotic environment and insect communities respond to 
these changes. This is supported by the fact that earlier studies 
have found regional genotype x environment interactions in the 
secondary metabolites of the same study saplings-^'. Yet, we do not 
know whether genotype x environment interactions in B. pendula 
resistance traits exist at a local scale and recent studies suggest that 
secondary metabolites are not the most important anti-herbivore 
defence of plants^ ^ On the other hand, spatial processes might 
affect local insect communities and create genotype x environment 
interactions. For example, in our experiment where genotypes of 
each block are arranged randomly, the effects of a particular geno- 
type could be partially masked by the effects of their conspecifics 
in some blocks if nearby genotypes are very dissimilar, i.e. there is 
associational resistance (see a review by Agrawal et al.^^') at the 
level of a genotype. Both of these processes may be affecting dif- 
ferent insect species differently. We found only one species pair 
that was correlated across genotypes in one of our study sites, 
which, together with earlier findings^^'^^ indicates that generalized 
defenses against multiple insect species are not likely in B. pendula 
(see Leimu and Koricheva^'O- Additionally, it may also be that local 
insect communities differ in their response regardless of spatial pro- 
cesses and without any change in the traits of B. pendula. 

The size of B. pendula trees is positively associated with their fitness, 
i.e. seed production^'^. It has been shown that herbivores can reduce 
the growth of B. pendula by up to 46% (Mikola et al. unpublished 



CN 

x 

03 
(/) 

Q 



O Natural herbivory 

□ Herbivore removal in previous year 



Stress 0.24 



23 




-2-1012 
nMDS axis 1 

Figure 4. Non-metric l\/IDS plots of insect assemblages of 27 
species colonizing 22 B. pendula genotypes that were either 
subjected to natural herbivory or protected from herbivory in 
the previous growing season in Kuikanniitty 2003. Each point is 
a centroid of six replicates. Numbers in the centre of the markers are 
genotype identification numbers. 



results, see also Prittinen et al.--, Silfver et al.-^') and increase 
seedling mortality considerably^^\ Thus, by imposing selection in 
various genetically variable resistance traits of B. pendula^^'^^'^\ 
herbivores may have high potential to drive the community evolu- 
tion in B. pendula. Indeed, we found that only one season of pro- 
tection from herbivory changed arthropod community variables 
(mean abundance and community composition) in five-year old 
field-grown B. pendula saplings. Total mean abundance, for exam- 
ple, was lower in saplings that were protected from herbivory in the 
previous growing season, which indicates that they may have had 
more resources to defend themselves against insects when herbi- 
vores were present again. Yet, the magnitude of these effects was 
smaller than the effects of local environmental (block) variation, 
and could explain only about 4% of the total variation in arthro- 
pod community structure. It is important to note, however, that in 
nature B. pendula seedlings typically establish in open patches, where 
high numbers of individuals compete heavily before self-thinning 
eliminates some of the seedlings. Surviving for these first years and 
consequently reaching maturity is crucial for an individual's fitness 
in this long-lived tree species. Earlier studies that have used open- 
pollinated progeny of the same genotypes, have shown that in such 
dense stands, even moderate levels of insect herbivory can change 
the genetic structure of B. pendula populations in the first year of 
establishment^-. This is reminiscent of recent studies, which have 
demonstrated that natural selection can favour different genotypes in 
the absence of herbivores rather than in their presence, and different 
genotypes in response to different herbivore species within only few 
generations of annual or biannual plants'' ^ ''^ (see also Hare^^). 
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To conclude, we have shown that the structure of insect herbivore 
communities can be significantly affected by intraspecific genetic 
variation when there is no mismatch in scale. However, genetic 
effects were modified by environmental variation on both a local 
and regional scale in one study year. Furthermore, insect herbivore 
damage in one growing season changed the community patterns 
of the following season, yet those effects were minimal compared 
to genetic and environmental factors. Our results suggest that both 
genetic and environmental factors are important determinants of the 
community structure of herbivorous insects. Together these mecha- 
nisms appear to maintain the high diversity of insects in B. pendula 
forest ecosystems. 

Data availability 

figshare: Community structure of insect herbivores on different gen- 
otypes of silver birch {Betula pendula), http://dx.doi.org/10.6084/ 
m9.figshare.91533256 
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This is an interesting study on the role of tree genotypes on herbivorous insect attack of willow trees, 
Betula pendula, using a comnnon garden approach in Finland. Twenty-two trees were randomly selected 
and cloned, although no details are given as to how the trees were cloned. Each tree is considered to be 
a different genotype. Six blocks with 4 saplings of each genotype were used from which one sapling was 
randomly selected to give 6 replicates per genotype for the study of insect attack. Insect sampling was 
done largely based on the damage done over the summer. I am surprised that it is possible to distinguish 
insects based on their damage as I would have guessed that the damage from leaf miners would have 
been similar, and damage from sucking insects difficult to find at all. I have the following questions and 
comments: 

1 . What is meant by structure of the insect community? This term is used in the abstract and 
throughout the paper but it is not defined. 

2. The term local environmental variation is mentioned, and it would be clearer if this was referred to 
as variation among blocks. 

3. Insects were removed from some saplings but details are lacking on how this was done, how 
frequently it was done, and how effective it was. 

4. "Additional Kuikanniitty 2003 data combined with Kuikanniitty 2003 raw data was used" I don't 
understand what this means. 

5. How would these results be interpreted by one who wanted to select for herbivore resistance 
among tree genotypes? Is the amount and consistency of the among genotype resistance 
sufficient to evolve over time? 

6. How are the results influenced by variation in insect abundance from year to year? 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard. 
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In this study, the authors selected and cloned 22 silver birch genotypes and planted them at two common 
garden sites to examine the effects of both genetic and environmental factors on the community structure 
of insect herbivores. By selecting genotypes from a spatially limited area (< 9 ha), and by performing the 
experiments at two spatially separated common garden sites, the authors were able to examine more 
precisely the main and interacting effects of genetic and environmental variation on the herbivore 
community. The authors also excluded herbivores on half of the host trees at each common garden site 
through application of an insecticide, which allowed them to examine the effect on the herbivore 
community in the following year. 

One limitation of the study is the use of only two common garden sites, which does limit the broader 
implications of the study. However, despite these limitations, the data still provide sound preliminary 
information on the importance of both genetic and environmental variation on the structure of the 
herbivore community. 

Although I understand and appreciate the challenges of documenting the herbivore community and the 
labor involved to do so, more comments are needed to address the lack of a more temporally robust 
sampling regime. Some species were sampled at a specific time given their respective seasonality, but 
most were sampled at the end of the growing season in the fall. In doing so, I suspect that authors would 
have missed spring and summer feeders, such as the winter moth for example, whose damage could 
have been difficult to ascertain and differentiate from other herbivores when herbivore damage was 
examined in the fall. This is not a fatal flaw, especially since the authors incorporated herbivore 
exclusions, which presumably would have also excluded spring and summer feeding herbivores. 
However, since one goal of the study was to examine the role of the herbivore community in affecting the 
community in the following year, it would be helpful to discuss the potential limitations of the fall sampling 
regime in documenting insects that feed earlier, and how the herbivore exclusion component was also 
(presumably) a mechanism to deal with sampling limitations. On a related note, it wasn't clear to me if 
potential temporal autocorrelation in the numbers and diversity of the herbivore community from one year 
to the next was appropriately addressed in the analysis; if so, I would suggest adding a statement in the 
materials and methods how this was addressed or if not, why it did not need to be addressed. 

Overall, this was a nicely designed experiment, the manuscript was very well written, and the data were 
well presented. This study adds to our knowledge of the role of genetic and environmental variation on the 
structure of the herbivore community on silver birch, and sets the stage for a number of interesting 
follow-up questions. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard. 
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